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Abstract 

The Color Glass Condensate (CGC) predicts the form of the nuclear 
wavefunction in QCD at very small x. Using this, we compute the wave- 
function for the collision of two nuclei, infinitesimally in the forward light 
cone. We show that the Wigner transformation of this wavefunction gen- 
erates rapidity dependent fluctuations around the boost invariant classical 
solution which describe the Glasma in the forward light cone. 
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1 Introduction 



The Color Glass Condensate (CGC) provides a description of the wavefunction 
of a hadron at very small values of x [1-5]. The CGC is a high density state of 
gluons which is controlled by a weak coupling, due to the high gluon density. 
Its properties are computable from first principles in QCD, at least in the limit 
of extremely high density. 

The CGC has also been applied to heavy ion collisions in order to generate 
the initial conditions for the evolution of matter in the forward light cone [6- 
10]. There are boost invariant solutions of the equations of motion supposedly 
appropriate for the high energy limit. The matter in the forward light cone, 
which we call the Glasma [10], initially has large longitudinal color electric and 
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magnetic fields, and as well as a large value of the Chern-Simons charge density. 
As time evolves, transverse color electric and magnetic fields appear as the 
remnants of the decaying longitudinal fields. This occurs as a consequence of 
the classical equations of motion. Eventually as the system becomes very dilute, 
these transverse color electric and magnetic fields may be treated as gluons, 
which form a quark-gluon plasma. The solution to the classical equations of 
motion is boost invariant and describes a system of expanding classical fields. 

It has recently been discovered that the boost invariant solution of the 
equations of motion are unstable with respect to rapidity dependent pertur- 
bations [11]. This instability has in fact close connections with the Weibel 
instabilities encountered in the physics of anisotropic plasmas [12-14], and it 
is speculated that such an instability may help the system created after heavy 
ion collisions reach a state of local equilibrium [15-21]. Whether or not such 
instabilities can grow to sufficient magnitude as to become as large as the clas- 
sical solution is the subject of current investigations. One piece of the puzzle 
which is not yet understood is the spectrum of the initial fluctuations. It is the 
purpose of this paper to derive an expression for the probability distribution of 
these fluctuations. 

If such fluctuations can grow to a magnitude comparable to the classical 
boost invariant field, then one has amplified initial quantum fluctuations to 
macroscopic chaotic turbulence. Such turbulence may ultimately be responsible 
for producing a thermalized Quark Gluon Plasma. Indeed, it proves difficult to 
thcrmalizc an expanding system of quarks and gluons by conventional multiple 
scattering formulae [22]. 

The evolution from the initial collision to the final state is shown in Fig. 1. 
Because of the similarities between the expansion of the universe in cosmol- 
ogy and the expansion of matter in heavy ion collisions, we call the latter the 
"little bang" . The initial singularity of cosmology is replaced by the singu- 
larity in the classical equations of motion associated with the collision. As in 
cosmology, where topological transitions associated with Chern-Simons charge 
may be responsible for generating baryon number, during the Glasma phase, 
there are also topological hclicity flip transitions [10,23]. It is during this time 
that instabilities develop, and there may be a Kolmogorov spectrum of density 
fluctuations generated [16,19,24]. This spectrum is similar to the spectrum of 
density fluctuations generated during inflationary cosmology. After inflation, 
the system reheats and thermalizes forming an electroweak plasma. In the little 
bang, thermalization might also be achieved after the Glasma expansion. Of 
course such an analogy between the little bang is not perfect, and it has not 
been established that the Glasma can in fact become thermalized. Neverthe- 
less, in the little bang, all of the physics is in principle understood from QCD 
in a small coupling regime, and one should be able ultimately to compute the 
evolution of the system. 

We begin this paper by a discussion of the evolution of an unstable sys- 
tem in quantum mechanics. We discuss a Gaussian initial wavefunction in the 
path-integral representation in which a description equivalent to the Schwinger- 
Keldysh formalism [25] naturally arises. We show that the distribution of initial 
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Figure 1: A light cone diagram illustrating the evolution of matter produced in 
heavy ion collisions 



fluctuations is given by the Wigner transform of the quantum mechanical wave- 
function [26]. These initial fluctuations are further evolved by solving classical 
equations of motion. We then proceed to show that the same procedure works 
in quantum field theory. We then compute the quantum field theoretical wave- 
function which describes the collision of two hadrons, at a time infinitesimally 
in the forward light-cone, after the collision. This should provide a practical 
algorithm for the computation of the evolution of matter produced in very high 
energy heavy ion collisions. 

The necessity of summing over the initial quantum fluctuations when solv- 
ing the classical Yang-Mills equations also appears when one goes beyond the 
classical approximation for the description of the fields in the forward light- 
cone. The formalism necessary for calculating loop corrections in a regime of 
strong sources and fields has been developed in Ref. [27] (and applied to the 
computation of the production of quark pairs - the simplest among the NLO 
contributions - in Ref. [28]). In this systematic approach, the instability of the 
boost invariant classical solution manifests itself as a divergence in the one-loop 
correction to the inclusive gluon spectrum, and it seems that the resummation 
of these divergent terms leads directly to an average over fluctuations of the 
classical initial conditions [29]. 

The description we advocate here has also some philosophical similarities 
with the work of Kharzeev et ai, who make an analogy to Hawking radia- 
tion [30]. We are unable to identify a Hawking temperature in our computation, 
but the idea that there is an initial spectrum of fluctuations generated at the 
initial singularity due to a quantum treatment of this singularity is similar. 
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2 The Case of Quantum Mechanics 
2.1 General Formulation 

We shall consider the time evolution of a generic quantum mechanical system 
from the initial time t 1 to the final time t F . If the system is as simple as an 
harmonic oscillator, the quantum fluctuations at the final time t F can be ex- 
pressed in terms of the initial quantum fluctuation at t t which evolve according 
to the classical equations of motion. In more complicated situations beyond 
the harmonic oscillator problem, the formulation we develop in what follows 
still holds unchanged, as long as we can make a Gaussian approximation in the 
weak coupling regime. 

Let us consider a simple example [31]. Supposing that we have a convex 
parabolic potential as sketched in Fig. 2 until the initial time tj , we know that 
the ground state has a Gaussian distribution of quantum fluctuations around 
the minimum of the potential. Let us now assume that the potential suddenly 
becomes inverted at tj : then the system is subject to an instability. If the 
problem was purely classical, the unstable state would not decay at all provided 
it lies exactly on the top of the potential curve. The decay into stable states 
is however unavoidable in quantum mechanics because of quantum fluctuations 
around the minimum. As long as higher order quantum fluctuations are negligi- 
ble, it should be acceptable to approximate the time evolution classically under 
such a potential causing the instability. As a result, we can anticipate that 
the convolution of the initial dispersion and the classical evolution provides the 
later distribution of quantum fluctuations at t F , and we shall verify this in the 
forthcoming discussion. We would like to emphasize here that our formulation 
is not limited to this sort of specific instability problem, but applicable to more 
generic problems. 



Figure 2: Schematic picture of the time evolution of quantum fluctuations. 
The ground state of a steady system for t < t I has a Gaussian distribution 
of quantum fluctuations. The distribution at a later time t F is given by the 
classical evolution from the initial distribution at t = t T . 

In the language of quantum physics the ground state wavefunction at t ; em- 
bodies the zero-point oscillation. The position and momentum fluctuate accord- 
ing to the uncertainty principle. The expectation value of physical observablcs 
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is defined as an average weighted by the wavefunction. Let us consider an ob- 
servable 0(x F ) which is a function of the position x F at the final time t F . [We 
specifically denote the position variable at t t and t F as x t and x F respectively 
in order to distinguish them from the position at other times.] By definition the 
expectation value at t F is given by 

(0)t F =jdx F \^(x F ,t F )\ 2 0(x F ) . (1) 

One can easily generalize the above expression to a situation involving observ- 
ables that depend not only on the position x F but also on the momentum p F . 
Using p F — —id/dx F and 0(x F , — id/dx F ) acting on ip(x F ,t F ), one can reach 
the same conclusion as in Eq. (10) in what follows, but for simplicity we only 
perform the derivation for the case of Eq. (1). 

The path-integral formulation naturally provides us the relation between the 
initial and final wavefunctions, 

if>{x p , t F ) = J dx, [Dx(t) Vp{t)] i>(x, , t, ) expj i jf F dt [px - H(p, x)] J , (2) 

where the integrals over the paths x(t) and p(t) extend from 1 1 to t F , with the 
boundary condition x(tj) = x 1 and x(t F ) = x F . The Hamiltonian is denoted 
by H(p,x). We will perform the path-integral approximately in order to find 
an analytically manageable expression. 

We expand the phase in the exponential around the stationary point, in order 
to reduce Eq. (2) to Gaussian functional integrals. The stationary condition 
leads to the classical path x = x c (t) and p = p c (t) determined by Hamilton's 
equations of motion, 

. _ 0ff(p°,s°) ._ dH{p°,x°) 

x ~ —dp~< ' p ~ (6 > 

with the initial conditions x c (t I ) — x z and p c {tj) = p t - We note that the final 
position x° (x r , p z ) = x c (t F ) is unique with x t and p z given, which means that 
we must consider the initial momentum p : as a function of x r and x F . Naturally, 
the quantum fluctuation around the classical path must be restricted at t 1 and 
t F to satisfy the boundary condition 5x(tj) = 5x(t F ) = 0. Up to the quadratic 
order in quantum fluctuations, therefore, we have 

ip(x F ,t F ) = J dx, ipix^tj) ex p|i^ dt[p c x c - H(p c ,x c )] 

J [D5x(t)VSp(t)] expj ij^ " dt SpSx ( 4 ) 



5x(t J )=5x(t F )=0 



op i ox z op ox 
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In the expansion of the argument of the exponential, the linear order terms in 
the fluctuations vanish because of the classical equations of motion. The first 
integral in the above expression is the classical part corresponding to the WKB 
approximation, where the x F dependence is implicit through the fact that x c (t) 
is the classical trajectory that starts at x 1 and ends at x F . The second integral 
represents the one-loop quantum corrections leading to the determinant of the 
propagator with the Dirichlet boundary condition. It is obvious from the above 
expression that, if H is the Hamiltonian of an harmonic oscillator (therefore 
consisting only of quadratic terms in x and p), e.g., H = ^p 2 + iw 2 i 2 , the WKB 
approximation is exact. It is because none of d 2 H/dp 2 ~l, d 2 H/dx 2 =ui 2 , and 
d 2 H/dpdx~0 depends on x I nor p T , and thus the integration with respect to 
fluctuations merely produces an irrelevant purely numerical factor. 

In general, the one-loop integral gives [det G^ 1 (x c )]~ 1 / 2 where G is the prop- 
agator of a fluctuation in the presence of the classical background x c . Usually, 
this determinant has no explicit p c dependence since usual quantum theories 
have at most quadratic terms with respect to the canonical momenta. These 
quantum corrections alter H{p c 1 x c ) to be H(p c ,x c ) — |lnG _1 (a; c ) and can be 
seen as an effective potential in the classical Hamiltonian. 

From now on we assume that the x c dependence in H{p c , x c ) is so strong that 
we can drop the x c dependence from the quantum corrections — ^ lnG~ 1 (x c ). 
This assumption should be checked case by case. Generally speaking, this is 
acceptable in instability problems like the one illustrated in Fig. 2. The quan- 
tum corrections modify the tree-level potential, which would eventually affect 
the equations of motion. But if the tree-level potential is steep enough, the al- 
teration due to quantum corrections should be negligible. Hence, our treatment 
is equivalent to assuming that the time evolution is dominantly governed by the 
tree-level potential. 

The probability distribution for the positions at final time is therefore cal- 
culated with the classical parts alone as 




The phase difference induced by the different initial conditions, x x and 2/ , can 
be simplified thanks to Hamilton's equations of motion (3) as 




(6) 
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From the middle to the last line only the surface term of the time integration 
remains. There is no finite contribution from the t — t F side because x F is 
common. It should be noted that p r (x T ,x F ) in the last line is a function which 
gives the value of the initial momentum p r as a function of the initial and 
final positions, i.e. p J (x / ,a; F ) = p I and p 1 {x' , x F ) =p' r Here, in order to get a 
convenient expression, we shall define P, as the initial momentum corresponding 
to the average of the two initial positions (and the same final position x F ) by 

P 1 ^p I {X n x F ) 1 (7) 

where X { = \{x I + x'j). Then we can write the last line of Eq. (6) as follows, 

r x i r\^ x i 
- \ d£ i Pi{x,,x F ) =- dA£ I p I (X I + Ax,,x F ) 

_ _J_ d 2 P I (X I ,x F ) 3 

- 1 1 24 8Xf 1 [> 

where we defined SXj = x r — x' t . Within the Gaussian approximation that we 
will employ later, it is enough for us to keep only the first term in the above 
expansion. 

The expectation value at the final time is now transformed in the WKB 
approximation from Eq. (1) into 

(0)t F = J dx.dx^dx; ^V,.0#S7>*/) e~ iP ' Sx ' 0{x F ) . (9) 

Here we can replace the x F integration by an integration over P t . Indeed, x F , 
Pj and Xj are related by the fact that the classical path that starts at the 
initial position Xj and with the initial momentum P T ends at the position x F , 
i.e. x F = x c F (X I ,P I ). Such a replacement is accompanied by the Jacobian 
Idx^/dPj] as a function of X n P n and t F . One can, in principle, obtain the 
Jacobian once one solves the classical path x c (x I ,p I ). This is, however, only 
a prefactor of the exponential for the phase of which we made the stationary- 
point approximation. It should be, therefore, consistent to drop the Jacobian 
prefactor off within the WKB approximation. 

After all, rewriting Xj = X t + \5x I and x' r = X T — ^5x I and integrating 
with respect to 8x I , we formally reach the following expression, 

(OK 

= J dX, dSx, dP^^X, - ISx^t^iX, + ISx^t^e-^i OixliX^P,)) 

= J dX I dP I WiX^P,) 0(a£ (X„P,)) , (10) 

where W(X I ,P I ) is the Wigner transform of the product of the two initial 
wavefunctions, defined as 

W{X n P I ) = JdSx, i)*{X I - \8x n t I )i}{X I + ±Sx I ,t I )e- iP i Sx ! . (11) 
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This is our formula and the purpose of our work is to calculate the Wigncr 
function analytically for the problem of heavy ion collisions, which should be 
a necessary input for numerical studies provided that the classical equations of 
motion are solved. 



2.2 Simple Example 

It would be instructive to demonstrate how the formula (10) works for an har- 
monic oscillator for which explicit calculations are feasible. Let us estimate 
the average of the observable O = x 2 F in a direct computation and by means 
of Eq. (10) and then compare the results. As we have already mentioned, the 
WKB approximation is exact in the case of the harmonic oscillator problem, 
which means that the two results should coincide. This is a trivial check of the 
validity of our formula (10). 

We shall characterize the initial state as a Gaussian [31]; 

^(x^t, = 0) = (irhy^ cxp(-x 2 /2b 2 ) . (12) 

The Green's function in the inverted harmonic oscillator with the Hamiltonian 
H = ip 2 — ^uo 2 x 2 , which is defined by 

(id t + d 2 + u 2 x 2 ) G(x, x'; t) = S(t)5(x - x') , 
lim G(x,x';t) = S(x - x') , (13) 

is known [31] to be 

G(x,x';t) = (27risinh^t)~ 5 cx p| 2 ~ ^ t {x 2 + x' 2 ) coshujt - 2xx' j, (14) 

from which we can directly calculate the wavefunctions at later time by 

ip(x F ,t F ) = J dxj G{x F ,x I ]t F )ip{x I ,Q) . (15) 

What we are calculating is the expectation value of the position dispersion 
when the potential is inverted at the initial time t } = 0. After some algebraic 
procedures we find, 

\^(x F ,t F )\ 2 = [^B 2 {t F )]-\^[-x 2 F /B 2 (t F )l (16) 

where B 2 (t F ) = b 2 cosh 2 ut F + b~ 2 smh 2 tut F . Then, it is straightforward to 
compute the average value of O = x 2 F , 

(0)t F = J dx F \i,{x F ,t F )\\l = \B 2 (t F ) . (17) 

This result indicates that the dispersion of the position grows exponentially with 
increasing time under the inverted potential just as sketched in Fig. 2. 
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Let us next consider the same problem according to the formula (10). The 
Wigner function obtained from Eq. (12) is 

WiX^P,) = i cxp(-X 2 /6 2 - 6 2 P 2 ) , (18) 

where we have adjusted the normalization factor, though it is only an unim- 
portant prcfactor, so that J dX T dP I W(X I , P T ) = 1. The classical path is 
fixed by Hamilton's equations of motion with the initial inputs x c (0) = X } and 
p c (0) — P n leading to the final point at t = t F as a function of X t and P n 

P 

xKX^Pj) = X I coshut F + ^smhcjt F . (19) 

After performing Gaussian integrations we can readily arrive at the expectation 
value. 

(0) tF = J dX.dP, ^(X J ,P J )4 2 (X J ,P J ) = l -B\t F ) . (20) 

We see that our formula perfectly reproduces the result from the direct evalua- 
tion in this simple case of the harmonic oscillator. The merit of our method is 
that we do not have to calculate the Green's function and the wavefunction at 
the final time directly. 



3 Generalization to Quantum Field Theories 

It is not difficult to extend the formula (10) described in Quantum Mechanics 
to the case of Quantum Field Theories. Let us consider in this section the case 
of a real scalar field theory. 



3.1 Introduction 

For a generic scalar field theory defined by a Hamiltonian density H[(f>, n] with 
field <f) and canonical momentum 7r, we expect that the counterpart of Eq. (10) 
is 

(0)t F = J [d&(*)] [cMsr)] Wfa,*,] O[#[0„7T,]] , (21) 

where (f) c F [(f) T , -k i ] is the solution of the classical equations of motion at time 
t F with the initial conditions (f>°(t I7 x) = <j)j(x) and 7r c (t J ,a;) = ir^x). This 
classical field is a solution of Hamilton's equations of motion, 

Here H denotes the integrated Hamiltonian H = J d 3 x Tt. The Wigner function 
is defined in the same way as in Quantum Mechanics by 

Wl^,^}^ J [dS^ix)] * , [^-i^ ) t J ]*[^ + ^ J1 i I ]e- i / d, " I ('WrW 

(23) 
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in terms of the initial wavefunction ^[(p,, t 7 ]. Similar expressions could be writ- 
ten for gauge fields after unphysical redundant degrees of freedom are removed 
by fixing the gauge. In the next subsection we will explain how one can construct 
the initial wavefunction W^,^] in practice. 



3.2 Wavefunction from the Annihilation Operator 

Let us consider the wavefunction and the associated Wigner function for the 
free real scalar field theory, which is the simplest extension to a Quantum Field 
Theory of the harmonic oscillator that we considered previously in Quantum 
Mechanics. In the Hamiltonian formulation, one can express the field operator 
(j)(x) and its canonical momentum operator n(x) in terms of the creation and 
annihilation operators as 

^ ) -/Svfc^ ci " + ^ c ^' (24) 

* (x) = - [ J w?^ k clkx ^ c " fe ' x ) ' (25) 

in the Schrocdinger representation in which operators are time independent. 
The dispersion relation for the free theory is u>k = \/ k 2 + m? . The canonical 
quantization condition is [4>(x),tt(x')] = iS^ (x — x') or equivalently [afc,aj,,] = 
(27r) 3 <5( 3 )(fc — k'). In momentum space, we have 

0(fc) = / d 3 * 4>(x) c- ik x = -^=(a h + a f _ k ) , (26) 
J V^k 

Tr(fe) = J d 3 x tt(x) c- ik x = -i^£(a k - ai fc ) , (27) 

for which the canonical commutation relation is deduced as [4>(k),n(k r )] = 
i(27r) 3 (5^ 3 ^ (k — k'). Let us consider an initial state which is an eigenstate of 
the operator </>(k), 

0(fc)*[&,t,] = 0,(fc) (28) 

From the commutator between 4>(k) and 7t(k'), we see that if(k) can be repre- 
sented as the differential operator with respect to the eigenvalue ^> 7 (— fe), that 
is, 

*(fc)*[^,t J ] = -i^-i^*[0 J ,t J ]. (29) 

Here the functional derivative in momentum space is understood in our con- 
vention as 5f(k)/6f(k') = (2n) 3 S^(k - k'). The ground state is the vacuum 
defined by the annihilation operator, i.e. 

o fc *[0 / ,t J ]=O. (30) 
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Thus, the ground state should follow from 



1 



'llOk (j), (fe) 



i,(-fe) 



*[(/> I ,t I ]=0. (31) 



It is easy to find the solution to the above equation being 



*[4>„tA 



d 3 fc 

(27T) 



3 ^^(-fe)^(fe) 



(32) 



The prefactor is the normalization for the Gaussian function, which is in fact 
unimportant in our approximation. The meaning of this wavefunction is that 
the ground state contains fluctuations with an energy w/c/2, that is, the zero- 
point energy. Then the Wigner function is 



W^^tTj] = iVexp 



r tfi 

J (27T 



d 3 fc_ 

) 3 



0U k 0, (-fc)& (fe) + — 7T, (-fc)7T 7 (fe) 



(33) 



with a normalization prefactor N independent of the fields and energy. One can 
clearly see that the result in the free real scalar field theory is just the prod- 
uct of the Wigner functions for an infinite assembly of independent harmonic 
oscillators (see Eq. (18)). 



3.3 Wavefunction from the Schrodinger Equation 

The strategy we used in the previous subsection, based on the canonical ap- 
proach, is obvious but not very suitable for more complicated problems like 
those encountered in gauge theories. We will develop here an alternate method 
in order to find the ground state wavefunction. 

In Quantum Mechanics, the wavefunction can usually be obtained as a so- 
lution of the Schrodinger equation when the Hamiltonian is provided. We will 
follow this procedure in the free real scalar field theory and will find a wave- 
function that is exactly the same as what we have obtained in Eq. (32). Let us 
begin with the Lagrangian density defining the free real scalar field theory, 

£ = - l(di4>)(dicP) - \m 2 <t> 2 ■ (34) 

The canonical momentum is it = 5L/6(dt4>) — dt<fr (where L is the integrated 
Lagrangian L = J d 3 x C) and the Hamiltonian density is by definition, 

H = nj> - C = ^tttt + h&iMdif) + \™ 2 cb 2 . (35) 

Of course, we could have started directly with this Hamiltonian density, but 
using the Lagrangian density as our starting point is more appropriate when we 
work with more general system of coordinates. This is so because the canonical 
momentum varies according to the choice of the time variable and so does the 
Hamiltonian density. 
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From the equal-time commutation relation, [<f>(t, x), n(t, x')] = 15^ (x — x'), 
we see that the action of the momentum operator tt{x) on an cigcnstatc of 
cj)(x) is identical to that of the derivative — iS/S(f)(x). The Schrodinger equation, 
idt^ = , can be expressed as a functional differential equation by means of 
this identification. The time dependence is actually separable by an ansatz \& = 
c -iEt-q, o w here *o is t independent. Then, the time independent Schrodinger 
equation arises as 



/ 



d 3 x 



+ -${x)(-d l d l + m 2 )<t>{x) 



2 5cj)(x) 2 ' 2 

The ground state solution to this equation is a Gaussian 



E^r 



(36) 



=N exp 



i J d 3 x (fij (x)y/— didi +m 2 4> I {x) 



d 3 k u k 
(2^T 



(37) 



which is identical to the wavefunction (32). Note that the energy eigenvalue E 
is given by 

f,, f d 3 fc Vfc 2 + m 2 

E =J ix Jw?^^- (38) 

which is nothing but the diverging zero-point energy. In principle, one could 
also write down the wavefunctions for excited states in terms of the Hermite 
polynomials. 



3.4 Wavefunction in terms of the r, rj coordinates 

For later convenience we shall see the same problem in the r, rj coordinates 
defined by 

re", x° - x 3 = re-" . (39) 



x° + x 3 



These coordinates are a natural choice for the purpose of describing the space- 
time geometry of a collisions between two high energy projectiles. Commonly r 
and rj are called the proper time and the rapidity respectively. The Lagrangian 
density in this coordinate system is 



re = ^(d T 0)(d T 0) - l-(d^)(d^) - • (0, 



-m 



2 i2 



(40) 



including the measure yjff^l = t. From this Lagrangian, we see that the 
canonical momentum is given by ir — rd T 4>. Then the Hamiltonian density is 



tTL = —tttt + — 
2t 2 



(£>„</>) (g>„0) + (0x0 • (9x0 + m 



2-2 



(41) 
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Now that we have the Hamiltonian density, we can write down the time inde- 
pendent Schrodinger equation as follows, 



J df] d 2 x_ L l^- 



1 



2t S^x^) 2 ' 2 
The solution of the above Schrodinger equation is then 



(42) 



*o[<£ 7 ] =A^exp 
=iVexp 



dk v d 2 fej_ 



(2tt) s 



'(Vr) 2 + 4+m 2 
^(-fc) V „ Mk) 



2t 



where we defined the Fourier transform, 



J d V d 2 X ± T <f>(T), x± ) e ~KknV+k±-x ± ) 



(43) 



(44) 



with fc = (k v ,k±). It is quite important to note that this wavefunction is 
t independent because the seeming r dependence is to be absorbed by the 
variable change rj — > rj/r in the first line and — > rfc^ in the second line. Of 
course, such rescaling changes the argument of the fields 7 (?7, x±) and </> 7 (fc) 
too. We will perform the integration over those fields at the end, however, 
and so the rescaling does not matter since possible t dependence in </> 7 (?7, x±) 
and <j>j(k) would be absorbed in this integration. This feature is essential for 
our construction of the wavefunction, for we made use of the time independent 
Schrodinger equation with an ansatz * = e~ iETl t>o where is t independent. 



4 Non-Abelian Gauge Theory 

We are now ready to proceed into the construction of the wavefunction in non- 
Abelian gauge theories. Before addressing the calculation of the wavefunction, 
we summarize our conventions, gauge choice, and basic equations. Then we 
will construct the wavefunction in the empty vacuum in the backward light 
cone region before the collision. We will connect the wavefunction from the 
infinitesimally backward light cone region to the infinitesimally forward light 
cone region with the boundary condition derived from the singularity in the 
Hamiltonian density. 
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4.1 Coordinates 

Let us first recall our conventions for the light cone coordinates and for the r, r\ 
coordinates 1 . The light cone coordinates are defined by 

x+ = -L(t + z) , x - = -L(t-z). (45) 

From there, we can readily obtain the metric tensor from the relation ds 2 = 

g^dx^dx", 

9+- = 9-+ = 1 , 9xx= g yy = -1 , (46) 

where the non-written components are all zero. Accordingly the light cone 
derivatives are to be understood as = d/dx T . 

The proper time r and the rapidity variable rj have already been defined in 
Eq. (39), and are related to the light-cone coordinates by x^ — (r / \/2)c ±n . The 
metric tensor for the r, r\ coordinate system is diagonal, with 

9tt = 1 , g m = -t 2 , g xx =g yy = -l . (47) 

Derivative with respect to r, r\ and x^ 1 can be related thanks to the following 
relation, 

x ± 

dx ± = — dT±x ± dr 1 . (48) 

T 

From the above identities, we can write the one-form gauge field 2 as, 

A T dr + A v dri = A~dx + + A + dx~ 

1 (49) 
= -(x+A~ + x-A+)dr + (x + A~ - x-A + )dr 1 . 

T 

Therefore, the components of the gauge field in the r, r\ coordinates are 

A r = A T = -(x+A- +x-A+) , A r , = -t 2 A n = x+ A- - x- A+ . (50) 
r 

In the same way, we can obtain the field strength as the two-form 3 F = dA — 
igA A A, that is, 

F Trj = d^-A^ — d ri A T — ig[A T , A^] , 

F Ti = d T A t - d t A T - ig[A T , A,] , 

F-qi = d^Ai — diA-q — \g\Aj), Aj\ , 

Fij = diAj - djAi - ig[A u Aj] , (51) 



1 Other definitions exist in the literature, that differ mostly in the way the light-cone coor- 
dinates are normalized. This changes the detailed expression of the metric tensor in light-cone 
coordinates. 

2 The transverse components have not been written because they are the same in all the 
coordinate systems we consider here. 

3 The symbol A denotes the anti-symmetric exterior product of forms. 



14 



where i refers to the perpendicular components (x and y) . We can raise or lower 
the indices using the metric tensor, e.g. F TV — g TT g m F T11 . One should keep in 
mind, however, the fact that F TV is not the same as d T A v — d v A T — ig[A T , A 71 ]. 
In order to avoid confusion, we shall work thoroughly in terms of quantities with 
lower indices, d T , d v , A T , and A n . The derivatives with the present conventions 
are 

8r = r=l, d v = -r>d^^, * = -* = ^- (52) 

4.2 Lagrangian and Hamiltonian 

In the r, r\ coordinates, the Lagrangian density is given by 

tC = —ixlcTg^F^F^] + r£ sourcc , (53) 



where we have included the Jacobian 4 \/\g^\ = r in the definition of the density, 
^source represents the singular source terms generated by the fast moving nuclei 
that collide at r = 0, 



-'source 



- Pl (x ± )S(x-)A- - p 2 (x ± )S(x+)A+ 



1 „ s 1 



-^-pi6(x ) + — P2&{x + ) A T - - -^—pi8{x ) '—P25{x + ) 

3j X > t + i" 



1 „, , 1 



A r) . 



(54) 



The source pi(x±_) represents the transverse color distribution of the nucleus 
moving to the positive z direction and the source p2(x±) is the same quantity 
for the nucleus moving in the negative z direction. These densities are random 
variables, and one should average observable quantities over the statistical en- 
semble of these distributions. All the considerations we make in the rest of this 
paper are for a given pair pi,p 2 taken in this ensemble. 

In the following discussion, we adopt A T = as our gauge fixing condition. 
This gauge has an important benefit in the problem we want to study. Indeed, 
one can see that it reduces to A + = if x+ = and to A~ = if x~ =0. 
This implies that the color current associated to the sources given in Eq. (54) 
is trivially conserved, because the color field produced by the collision is such 
that the terms that may have induced a precession of the current ([A + , J~] 
or [A~ , J + ]) arc identically zero. In practice, this means that we can safely 
consider the color current as given once for all in terms of pi and p 2 , and forget 
about the current conservation relation. 

In fact, naively, because A~ — > when x~ — > and A + — ► when x + — ► 0, 
one may have the impression that the whole source term in the Lagrangian 
density has no effect at all. However, this conclusion is incorrect, because one 
has to differentiate with respect to the fields in order to obtain the equations of 
motion, the latter have a contribution from the source terms. Moreover, in the 
r, 7] coordinates, one has 

A v - C(r 2 ) -» (55) 



4 One has d 4 x = rdrdrid 2 x±. 
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when t — > 0, as one can see from the definition (50). In terms of the field 
component in these coordinates, the Lagrangian density in the A T = gauge 
reads 



tC 



-tr 



7 (d T A T ,) 2 -2(d T A i f 



2 ^vi + FijFi; 



source j 



(56) 



and the canonical momenta are given by 

d(rC) 



El d{d T A t ) 



= rd T Ai = -rF h 



1 



1 



(57) 



d{d T A ri) 

As a consequence, the Hamiltonian density is 

tH =2tr[(d T A l )E t + (d T A v )E' n ] - tC 
=tr 



7]T 



-E l E l + tETIETI + -F m F m + 

T T 2 



-^—pi5(x ) —p 2 5(x + ) 



(58) 



In the case of gauge theories, the time evolution is not uniquely generated by 
the Hamiltonian due to the redundancy of the gauge degrees of freedom. It 
should become unique when the Gauss law constraint, which is deduced from 
the coefficient of A T in the Hamiltonian density (before the gauge condition is 
enforced), 



DiE* + D, q E' n - — 



-^-piS(x ) + —p 2 5(x+] 



o, 



(59) 



is satisfied by the physical state. The last term of the Gauss law looks like it may 
be zero (thanks to t 2 = 2x + x~), but we should keep it because it contributes 
to the singularities of the Hamiltonian when divided it by r, from which we can 
determine the boundary conditions for the classical fields as shown in the next 
subsection. 



4.3 Classical background fields 

From the Hamiltonian, H = f drj d 2 x±H., one gets the classical equations of 
motion for the classical background fields A and momenta £ , which should also 
satisfy the Gauss law (59), as [11] 



d T E l 



<9 T £" = 



S(tH) 
SAi 

S(tH) 
SA n 



-Djj-jjj 



: pi5(x~ 



:P26(x+] 



(60) 
(61) 
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T evolution 




const. T| 

Figure 3: The source singularities are shifted by an infinitesimal amount e in the 
x + and x~ directions respectively. The r integration from 0~ to is taken 
along a line of constant r\. The backward light cone region corresponds to a 
negative r. 



with the background covariant derivative T> = d — igA, and the electric fields 
£ 1 = rd T Ai and £ v = T~ 1 d T A n . These equations must be solved with retarded 
boundary conditions, in which the fields and momenta are all vanishing when 
t — ► — oo. The solution of these equations of motion is known analytically in the 
space-like regions below the forward light-cone. Inside the forward light-cone, 
the solution must be obtained numerically, with an initial condition at r = + 
which can be obtained analytically. In order to determine this initial condition, 
one must integrate the equations of motion in an infinitesimal region just above 
the forward light-cone. 

In terms of the sources p\ and p%, these initial conditions are known to be : 

Ai = a] + a\ , (62) 
A v = , (63) 
£' = rd T Ai = , (64) 

= ~d T A n = ig[a\,o%], (65) 

where a\ and are known functions that depend only on the perpendicular 
coordinates x±, and satisfy dia\ — p\ and diCtf = p2- (The repeated index 
i indicates a sum over the transverse indices x and y.) In the rest of this 
subsection, we rederive these initial conditions in the r, rj coordinates. This is 
instructive as it illustrates some of the issues encountered when using the r, 77 
coordinates, and it will pave the way for the forthcoming construction of the 
wavefunction. 
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The main difficulty when dealing with the r, t] coordinates is that the (physi- 
cal) singularity of the sources is located at the same position as the (unphysical) 
coordinate singularity 5 at t = 0. It should be noted that, while the backward 
light cone region is parameterized by a negative r, the difficulty stems from the 
impossibility to describe the space-like regions by a real r. We will introduce a 
regulator e, that we use to slightly shift the location of the sources. This helps 
to separate the physical singularities from the unphysical ones, as illustrated in 
Fig. 3. Then the space-time is divided into four regions, 



Region I 
Region II 
Region III 
Region IV 



x + < e, x~ < e, 
< x + < e, x~ > e, 
< x~ < e, x + > e, 
x + > e, x~ > e, 



separated by the source singularities. Let us consider the r evolution along a 
line of constant 77, e.g. in the negative 77 region as illustrated in Fig. 3. Then 
the r evolution first goes through a singularity at x~ = e from the region I to 
the region II, and then through a second singularity at x + = e that leads to the 
region IV. 

Equation (63) is trivial from Eq. (55). By inserting the canonical momenta 
into the Gauss law divided by r, we can get the following boundary condition, 



C 



—8 8 A — - 

r 2 U V U T-^V 2 



■\p 1 8{x -e) + — p 2 S(x + - 
x + x 



, (66) 



where we perform this integration (66) with respect to r in the small range 
from CP to ( at fixed -q. The upper bound r = ( should be chosen such that 
( > v / 2ee ?? and £ < V2ee~ v , so that the integration ends in the region II 6 . 
Since only the singularity at x~ — e is present on the integration path, we can 
assume that the background fields in the region II are proportional to the step 
function 9(x~ — e). Note that such a step function gives a delta function when 
differentiated with respect to r or 77. Let us therefore substitute the following 
forms, 



At = a\6{x- - e) , A v = T 2 (3 1 e{ x - - e) 



(67) 



into Eq. (66). The only terms that survive in this integral in the limit e — > + 
are those containing a delta function. It should be noted that the Ai in the 
covariant derivative T>i drops out because it has the same color structure as 
d T Ai, as long as only the singular part is concerned (obtained when d T acts on 
the step function rather than on a}). After all, we get 



dta] - 2/3 1 - Pl = 



(68) 



5 The singularity of this system of coordinates at r = can be seen in the fact that the 
determinant of the metric vanishes at r = 0. 

6 Thanks to our shift of the sources by the amount e, this path is entirely contained in the 
time-like region, and can be parameterized by the r, rj coordinates. 
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The equation of motion (60) is then satisfied almost trivially, and the other 
equation of motion (61) leads to the very same boundary condition, and there- 
fore we cannot determine uniquely aj and j3\ . This degeneracy is related to the 
residual gauge freedom; if we allow a nonzero 1 proportional to pi, then we 
can eliminate it by a gauge rotation 7 applied to A n without affecting the gauge 
condition A T = (except for terms of order r — > 0) . We shall choose fi 1 = and 
therefore diOt\ = p\. We can do the same for the region III, where we obtain 
the background Ai = a 2 9(x + — e) with diaf = p 2 . 

Then, we shall push the time integration towards the region IV. There is 
another singularity at x + = e for rj < and at x~ = e for rj > 0, which can be 
accounted for by using the following parameterization, 

A l = a\6{x~ -e)6{e-x + ) + a 2 i 6{x + -e)6{e-x~) + a\6{x~ -e)6{x + -e) , (69) 

where a\ and a 2 are already fixed. If we differentiate twice with respect to r or 
T] the products of theta functions in this expression, we obtain terms with two 
delta functions that are singular at t = rj = (in the limit e — > + ). From the 
equation of motion (60), we therefore get the condition, 

f o+ f C i 

/ dry / dr-^A = 0, (70) 
Jo- Jo- T 

where this time we must chose £ > y/2ee ±v so that the endpoint of the r 
integration is in the region IV. This relation leads to 

aj = a] + of . (71) 

Let us finally perform the integration (66) with the ansatz 

A v = T 2 f3 3 6(x + - e)6{x- - e) . (72) 

In this case, it is important to note that d T Ai can have a different color structure 
than that of Ai . Therefore, we should keep the covariant derivative T>i and treat 
carefully the commutators such as [Ai, d T Ai] and [Ai, d^Ai). The Gauss law (66) 
eventually leads to the boundary condition 

d l a\ + d l a\ - ig[a\ , a 2 ] + 2/3 3 - Pl - p 2 = , (73) 

and the equation of motion (61) gives again the same condition, which results 
in 

s = a?]. (74) 

Hence, one can reproduce the classical background fields (62)-(65) in the r, r\ 
coordinates by slightly shifting the source singularities. We should emphasize 
here that, even though we have £ l = 0, we must keep terms such as d T £ l and 
t~ 1 £ 1 which are non-vanishing when r — > 0. Since A v , in contrast, goes to zero 
like 0{t 2 ), we can drop d T A ri = and t~ 1 A, 1 = in the limit r — > 0. 

7 This situation is actually analogous to what happens when solving the equations of motion 
in the light-cone coordinates in the A + = gauge, in the presence of the source singularity 
5 fl +pi5(a; _ ). In this case, A~ can be gauged away by an x + independent gauge rotation. 
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4.4 Wavef unction in the backward light cone region 



We now proceed with the derivation of the wavefunction in the r, r\ coordinates 
in the non-Abelian gauge theory. Let us first consider the backward light cone 
region in the weak coupling regime. We see that A v and E n are not a dynamical 
variables on the light cone in the A T = gauge. Therefore, the wavefunction 
does not involve A v at all and the Gauss law provides E v as a function of the 
dynamical variables. In other words, E n is a solution of the Gauss law with 



A, q = substituted, 



Or, 



(75) 



where we should remark that the fields and the momenta in this identity are the 
full ones, including both the classical backgrounds and the quantum fluctuations. 
In the following, we denote the quantum fluctuations with lowercase letters, 



di — Ai Ai 



®r] — A 7 ] Ajj , 



e { = E l -£\ e v = E 71 - . 



(76) 



We then expand the Hamiltonian to quadratic order in the fluctuations, and 
then the time-independent Schrodinger equation can be written as 



J dr/ d 2 £Cj_tr| 



1 <5 2 S 2 1 



Sal 



+ 



(Vidj - Vjdi) 2 - 2\gFij [a,, aj \ \V- [A] = [A] , 



(77) 



with an energy eigenvalue E. When writing this Schrodinger equation, we have 
replaced the canonical momenta by differential operators, thanks to the canon- 
ical commutation relations at equal t, 

[aUv,x±),e j Hv\x'±)]=iS ab S ij S(ri-ri')S {2 Hx±-x'±) , (78) 
[a a v ( V ,x ± ),e" b ( V ',x' ± )] = iS ab 5( v - V ')S^ (x ± -x' ± ) . (79) 

It should be mentioned here that the right hand side does not involve r explicitly 
because we have defined the canonical momenta e l and e n from tH, including 
the measure. 

Changing temporarily to the variables = a n /r (and d v = rd^), and using 
the notation a 2 = (a i; a^), we can write the Schrodinger equation in the following 
compact form, 



/ 



d^d x± tr 



t 5a 2 , 



+ ra : ■ G } ■ a j 



*_L4] = EV_[A] 



where we have defined the background propagator, 

{G-}) ab = —5 IJ (D 2 ) ab + {V 1 V J ) ab - 2gf acb F c I . 



(80) 



(81) 
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with T> 2 = d 2 + T>iDi. The formal solution to this Schrodinger equation is 
not difficult to find. The ground state solution to Eq. (77) or Eq. (80) can be 
expressed in terms of the square- root of the inverse propagator G^ 1 as 



4»_ 



N exp 



J drj^x^ t tr (a, <Jg~J a,^ 



(82) 



where we have chosen only the normalizable solution. Since there is no back- 
ground field in the backward light cone region, the square-root of the inverse 
propagator is easy to calculate. It is simply proportional to the transverse pro- 
jection operator d IT — djdj/d 2 . Noting that = r -1 ^ — > 0, we see that the 
wavefunction is simply 



[A] = N exp< - / dr] d 2 x± tr 



(li T 




did 



(83) 

Except when d v is smaller than tQ s , the d v /r term dominates over the trans- 
verse derivative. For the analysis of the instability, however, we may have to 
retain the transverse momentum here, because the unstable rj dependent modes 
may be very soft [11]. Also we note that, although we have used the notation Oj, 
it is identical to the full fields Ai in the absence of backgrounds in the backward 
light cone region. 



4.5 Singularity in the Hamiltonian 

In order to clarify the boundary condition at the source singularity, we first need 
the Hamiltonian in the presence of background fields, and more specifically its 
singularity. Note that the singular part of the Hamiltonian should commute 
with the Gauss law because of gauge invariance. 

The explicit form of the Hamiltonian in the presence of the background reads 



tH = tr 



-eV + -£ i e i + -£ l E l + re^e" + 2 T £^e^ + t£ t '£^ 

T T T 



+ - (dndi -Vicir,) 2 \g (9„A)[a„' a *] + -(^A)^^ -X>jO„) + -{dnAif 

T T T T 

+ 2" { v i a 3 - Vjdi) 2 - rig T io [a,, a,-] + t T lj (d l a j - djcn) + 



+ r\^-p 1 S(x ) —p 2 5(x + )\i 

lx + x J 



(84) 



up to quadratic order in the fluctuations. The last term is the current density 
due to the sources, which can also be expressed as J n = —(2/r)did v Ai. Gauge 
invariance requires the covariant conservation of the current, i.e. D^J 71 = 0, 
which is actually satisfied thanks to Eq. (70), since A v — > at r — > . 



Thus, in the limit a v — ► 0, the singular terms in the Hamiltonian are 



tr 



2{d T Ai)e l + -(d v Ai)(d v a,i) + -{d^f 



(85) 
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Here we have used Eq. (57) in order to rearrange the first term. The delta- 
function singularities originate from d v and d T acting on the step function in 
A{. The second and third terms can be in fact dropped thanks to Eq. (70), after 
an integration by parts. Therefore, the only remaining singular term is after all, 

rising = 2tr[(d T A)e l ] , (86) 

on the light cone where a v is vanishing. 



4.6 In the forward light cone region 

By integrating the Schrodingcr equation with respect to r from 0~ to £ with 
the singularity in the Hamiltonian we have just identified in Eq. (86), we can 
derive a relation between the wavefunctions in the infinitesimally backward and 
forward light cone regions. 

In principle, we can follow the same procedure as the one used for the deter- 
mination of the classical backgrounds. By shifting the singularities as shown in 
Fig. 3, we can access the space-like regions while using the r, ij coordinates. In 
the presence of the regulator e, we should keep A v terms which are of order e 2 , 
and the r integration from I towards IV will pick up singular terms proportional 
to A v . When the wavefunction contains linear terms in A v in the exponential, 
the A v integration in the Wigner transformation results in the delta-function 
constraint imposed on the conjugate variable E v . At the end in the limit of 
e — > 0, we must take simultaneously the limit A v — > 0, so that the A n terms in 
the wavefunction disappear. Nevertheless the constraint on E n preserves some 
information about the singular terms. This procedure is a bit tedious to imple- 
ment, but the final answer is quite simple: the Gauss law always holds, from 
which one can obtain E n as a function of the other fields. Therefore, the con- 
straint on E v after the A v integration simply coincides with the Gauss law (59). 
This means that we can obtain the final answer by the following shortcut: drop 
all the singular terms involving A v and add by hand to the Wigner function a 
delta-function constraint fixing E v from the Gauss law. 

Since we can set A v — even with finite e, the singular Hamiltonian (86) is 
sufficient for our purpose to connect the wavefunctions between the two light- 
cones. The Schrodinger equation leads to the boundary condition from r = CT 

tO T = C -> 0+, 



i(*+-*_)= dr J d V d 2 Xl _ rW sing )i (*+ + *_) . 



(87) 



Neglecting higher order terms in the fluctuations, we can solve the above equa- 
tion as 

^+ ~ cxp — ij dr J drjd 2 x±TH s 



*- , (88) 



where TH s i ng is given by Eq. (86). The identification of e l as the differential 
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operator —iS/SAi enables us to write 



Jo- 



f_dTTH sins = 2tr[A(CK] = -^i(C)^ • 



(89) 



The exponential of this differential operator is nothing but the translation op- 
erator of Ai by an amount — A\ ; i.e. Ai — > Ai — Ai. Therefore, we conclude 
that V + [A] =^_[A — A]. By taking the limit e — > 0+ while keeping the final 
time £ very small but fixed, the surface of constant time r = C lies entirely in 
the region IV, where the classical transverse field is Ai = aj + a 2 . Thus, we 
finally obtain the following wavefunction in the forward light-cone, 



with the wavefunction defined in Eq. (83). 

Although the wavefunction takes the same functional form, it is the Gauss 
law which is actually affected nontrivially by the source singularity, and the 
fluctuations in e v are uniquely fixed from the background fields (62) and (65), 
and the fluctuations a, and e 4 . Namely, from Eq. (75), up to term of quadratic 
order in fluctuations, we have 



where we assumed that the fluctuation e l does not contain a singularity of the 
type ~ Td T 9(x ± — e) (unlike the classical transverse electric field E % - The 
r/ integration of T>iE % leads to £ v because of this singularity.) The origin of 
the presence of an indefinite integral constant is that the Gauss law cannot 
constrain the r\ zero-mode of e n . In such a case, the a v integration in the Wigncr 
transform would impose the zero-mode to vanish. Hence, e^J should be fixed by 
the condition f^^drje 71 = 0. 

4.7 Wigner transformation 

Finally, we can compute the Wigncr function from the wavefunction. The con- 
ditions of a v — > and the Gauss law constraint are, as we have noted before, 
expressed by delta-functions added by hand. We can then write the Wigner 



*+L4] = * _[A - {a 1 + a 2 )] , 



(90) 




(91) 
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function as follows 



W[a,e] 

= N f [dSa] **[a-|<ya] *+[a + ±<5a] C - [ f d vd 2 x± 2tr[ e i 5a i ] 



iVexpi ~y dr;d 2 a;_ L 2tr ai t ^ - (d v / t) 2 - d]_ (<% 



did j 



V lJ + (dJrY) 



r\/-(V) 2 -S± 



x <5 [a,,] <5 



e^-el+ drj V { e 



(92) 



where Sij + didj / (d v / t) 2 is the inverse of J y - & 0,- / [ (0„ It) 2 + d\] . This Wi gner 
function allows us to compute the average of any observable 0[A] at later time 
r > with initial quantum fluctuations, as follows 8 , 



(0) T = J [da, da„ de l de"] W[a, e] O [^[A + Oi, a v , e\ £ " + e r '; r]] 



(93) 



Here we symbolically denote by A[Ai + o», a,,, e\ f + e' 7 ; r] the classical gauge 
field obtained at time r by solving the classical equations of motion (60) and 
(61) with initial conditions A% + cii, a v , e\ and E 71 + e v at r = + , where Ai 
and are the usual background fields (62) and (65) respectively. 

The Gaussian Wigner function describes the dispersion of fluctuations around 
the classical initial condition. The distribution turns out not to be of white- noise 
type but to be characterized by the correlation function, 



(a i (r],x±)a j (r]',x' ± )) 



J [da,i da v de 4 de 1 '] W[a, e] ai(rj,x±) a,j(r}' ,x' ± ) 



:(% + 



did* 



[da, da v de 1 de n ] W[a,e] 



r^J '-(V) a - dV ' (9 " /T)2 
The fluctuation of canonical momenta is characterized inversely by 
(e i (??,x_ L )e j (?7',a;' ± )) 



(94) 



= T^-(d v /r) 2 -d\(5 l0 



did. 



(d v /r) 2 + d[ 



(95) 



8 A result very similar to ours was already known in the case of a scalar field subject to a 
parametric resonance [32] ; namely that calculating the expectation value of some combination 
of fields can be done with a solution of the classical equation of motion, provided one averages 
with a Gaussian weight over fluctuations of the initial condition for this classical solution. See 
also Ref. [26]. 
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The formula (93) and the above two relations, Eqs. (94) and (95), are our 
central result. Regarding the rj components, a v is vanishing while e n is fixed 
by the Gauss law. In the limit of r — > the distribution becomes simpler, 
and depends solely on the derivative with respect to rapidity, d^. For practical 
uses, however, it is not possible to start the numerical simulation exactly at 
r = + , and the transverse derivatives d± will then play the role of a regulator 
for infrared modes when d v /r is smaller than d±. 

5 Conclusions 

In this paper, we have discussed an approximation scheme for calculating the 
quantum expectation value of an observable that depends on field configurations 
at late times. In the WKB approximation, we have obtained a formula giving 
this expectation value from the classical equations of motion, with an initial 
condition that includes fluctuations. The spectrum of these initial fluctuations is 
specified by the Wigner transform of the wavefunction, the latter being obtained 
as a ground state solution of the Schrodingcr equation. Our formula is generic 
enough that we can expect a variety of applications, for instance to the study 
of an unstable system, or to that of the spinodal decomposition associated with 
a phase transition. 

An application of particular interest to heavy ion collisions is to calculate 
the Wigner function of these fluctuations immediately after the initial impact of 
the two nuclei, in the weak coupling regime. These initial quantum fluctuations 
provide the seeds that trigger the instability of the rapidity dependent modes. 
The precise determination of initial fluctuations or seeds is of crucial importance 
for the estimate of the growth rate of the unstable modes. 

By identifying the physically relevant gauge fluctuations and the singular- 
ity in the Hamiltonian density, we have derived the discontinuity between the 
wavefunctions in the backward and forward light-cones. It turns out that the 
initial singularity results only in a shift in the transverse gauge fields, equal to 
the classical background field. This is because, on the light-cone, the physical 
degrees of freedom are gauge fields polarized in the perpendicular plane and the 
longitudinal mode along the collision axis has to vanish. 

This wavefunction and the resultant Wigner function describe the spectrum 
of initial fluctuations. There are two noteworthy features in our results. 

Firstly, both e 1 and dj fluctuate, with 2-point correlations that are inversely 
related to one another. This means, if the boost invariance is largely violated by 
e l , then ai fluctuations are suppressed, and vice versa. Thus, the prescription 
adopted in Ref. [11], in which only the e 1 fluctuations are taken into account, 
needs to be extended in order to include also the fluctuations of ai . It should be 
important to describe these fluctuations correctly because the ai distribution 
rather than e l distribution has large correlation for infrared modes having a 
small rapidity dependence. 

Secondly, although it may be less important, the off-diagonal correlation 
between x and y components is not vanishing. Again, this might be relevant for 
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infrared modes. This feature is also missing in Ref. [11]. 

It would be interesting to see how these modifications of the fluctuation 
spectrum affect the growth rate of the instability. Although the spectrum of 
the initial fluctuations could be obtained analytically, this question can only 
be investigated numerically since it requires to solve the classical Yang-Mills 
equations in the forward light-cone. 
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